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Abstract 

We evaluate the kinematic distributions in phase space of 4-parton final-state subprocesses pro- 
duced by double parton scattering, and we contrast these with the final-state distributions that 
originate from conventional single parton scattering. Our goal is to establish the distinct topolo- 
gies of events that arise from these two sources and to provide a methodology for experimental 
determination of the relative magnitude of the double parton and single parton contributions at 
Large Hadron Collider energies. We examine two cases in detail, the b b jet jet and the 4 jet final 
states. After full parton-level simulations, we identify a few variables that separate the two con- 
tributions remarkably well, and we suggest their use experimentally for an empirical measurement 
of the relative cross section. We show that the double parton contribution falls off significantly 
more rapidly with the transverse momentum p£ of the leading jet, but, up to issues of the relative 
normalization, may be dominant at modest values of pfi ■ 
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I. INTRODUCTION 



Double parton scattering (DPS) means that two short-distance subprocesses occur in a 
given hadronic interaction, with two initial partons being active from each of the incident 
protons in a collision at the Large Hadron Collider (LHC). The concept is shown for illustra- 
tive purposes in Fig. [TJ and it may be contrasted with conventional single parton scattering 
(SPS) in which one short- distance subprocess occurs, with one parton active from each initial 
hadron. Since the probability of single parton scattering is itself small, it is often expected 
that the chances are considerably suppressed for two or more short-distance interactions in a 
given collision. However, expectations such as these bear quantitative re-examination at the 
LHC where the high overall center-of-mass energy provides access to very small values of the 
fractional momentum x carried by partons, a region in which parton densities grow rapidly. 
A large contribution from double parton scattering could result in a larger than otherwise 
anticipated rate for multi-jet production and produce significant backgrounds in searches for 
signals of new phenomena. The high energy of the LHC also provides an increased dynamic 
range of available phase space for detailed investigations of DPS. 
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listory theoretically 



21 



22], and there is evidence for 



their presence in collider data from the CERN Intersecting Storage Rings [23| and Fermilab 
Tevatron 
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251 ] . A significantly greater role for double-parton processes may be expected 



at the LHC where higher luminosities are anticipated along with the higher collision energies. 
Of substantial importance is to know empirically how large the double parton contribution 
may be and its dependence on relevant kinematic variables. 




FIG. 1: Sketch of a double-parton process in which the active partons are i and k from one proton 
and j and I from the second proton. The two hard scattering subprocess are A{i j — > a b) and 
B(kl^c d). 
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Our aim is to calculate characteristic final states at LHC energies in which it may be 
straightforward to discern a double parton signal. We show in this paper that double 
parton scattering produces an enhancement of events in regions of phase space in which the 
"background" from single parton scattering is relatively small. If such enhancements are 
observed experimentally with the kinematic dependence we predict, then we will have a 
direct empirical means to measure the size of the double parton contribution. In addition 
to its role in general LHC phenomenology this measurement will have an impact on the 
development of partonic models of hadrons, since the effective cross section for double parton 
scattering measures the size in impact parameter space of the incident hadron's partonic hard 
core. 

From the perspective of sensible rates and experimental tagging, a good process to exam- 
ine should be the 4 parton final state in which there are 2 hadronic jets plus a b quark and 
a b antiquark, viz. b b j\ jz- If the final state arises from double parton scattering, then it is 
plausible that one subprocess produces the b b system and another subprocess produces the 
two jets. There are, of course, many single parton scattering (2 to 4 parton) subprocesses 
that can result in the b b j\ ]2 final state, and we look for kinematic distributions that 
show notable separations of the two contributions. As we show, the correlations in the final 
state are predicted to be quite different between the double parton and the single parton 
subprocesses. For example, the plane in which the b b pair resides is uncorrected with the 
ji 32 plane in double parton scattering, but not in the single parton case. 

The state-of-the-art of calculations of single parton scattering is well developed whereas 
the phenomenology of double parton scattering is as yet much less advanced. In the remain- 
der of this Introduction, we first describe the approach we adopt for the calculation of double 
parton scattering, specializing to the proton-proton situation of the LHC. Then we outline 
the paper and summarize our main results. Our calculations are done at leading-order in 
perturbative QCD, adequate for the points we are trying to make. 

Making the usual factorization assumption, we express the single-parton hard-scattering 
differential cross section for p p — > a b X as 



i,3 

Indices i and j run over the different parton species in each of the incident protons. The 
parton-level subprocess cross sections da^j^ a ^(xi, x[, /i) are functions of the fractional par- 




(1) 
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tonic longitudinal momenta x± and x[ from each of the incident hadrons and of the partonic 
factorization/renormalization scale /x. The parton distribution functions fp(xi,/i) express 
the probability that parton % is found with fractional longitudinal momentum x\ at scale 
H in the proton; they are integrated over the intrinsic transverse momentum (equivalently, 
impact parameter) carried by the parton in the parent hadron. 

A formal theoretical treatment of double parton scattering would begin with a discussion 
of the hadronic matrix element of four field operators and an explicit operator definition 
of two-parton correlation functions. This procedure would lead to a decomposition of the 
hadronic matrix element into non-perturbative two-parton distribution functions and the 
corresponding hard partonic cross sections for a(ijkl — > abed). An operator definition of 



26] where the two-parton correlation 



two-parton correlation functions may be found in Ref. 
function is reduced to a product of single parton distributions. An explicit operator definition 
of two-parton distributions with different values of the two fractional momenta Xi and x<i 



is presented in Ref. 22J, along with a model for the two-parton distributions in terms of 



normal parton distributions. In this paper, we fol 



lines similar to Refs. 
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In a double parton process, partons % and k are both active in a given incident proton. We 
require the joint probability that parton k carries fractional momentum X2, given that parton 
i carries fractional momentum x\. In general, this joint probability H t,k (xi, X2, fiA, /Ab) 
should also depend on the intrinsic transverse momenta fc*r,j and kx,k of the two partons 
(or, equivalently, their impact parameters). The hard scales ^a and fiB are characteristic of 
the two hard subprocesses in which partons i and k participate. In the sections below, we 
discuss the choice we make of the hard-scale and do not explore in this paper theoretical 
uncertainties associated with higher-order perturbative contributions. 

In contrast to single parton distributions functions fl,{x\, /i) for which global analyses have 
produced detailed information, very little is known phenomenologically about the magnitude 
and functional dependences of joint probabilities H l,k (x\, X2, fiA, Hb)- A common assump- 
tion made in estimates of double parton rates is to ignore possibly strong correlations in 
longitudinal momentum and to use the approximation 



H l ' k (x 1 ,X 2 ,IJ,A^B) = rJx u fi A )fp(x2,fJ lB )- 



(2) 



For reasons of energy-momentum conservation, if not dynamics, the simple factorized form of 
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Eq. (j2J) cannot be true for all values of the fractional momenta x. The values of £2 available 
to the second interaction are always limited by the values of x\ in the initial interaction since 
X\ + £2 < 1- The approximation certainly fails even at the kinematic level if both partons 
carry a substantial fraction of the momentum of the parent hadron. However, it may be 
adequate for applications in which the values of X\ and £2 are small. We remark that the 
momentum integral 

E/W^*,,*,,,,,,,.)^-!, (3) 

i,k J 

as long as we can run the upper limits of the x\ and £2 integrations to 1, independently. 
The large phase space at the LHC may make it possible to explore dynamic correlations 
that break Eq. 0. 

In Fig. [21 for the region of phase space of interest to us, we show the contributions to 
the bbjj cross section as a function of x from both DPS and SPS, after minimal acceptance 
cuts are imposed (Sec. [TTJ) . The center-of-mass energy is ^Js = 10 TeV. It is evident that 
the majority of DPS events are associated with low x values, in essence never exceeding 
0.2. The momentum carried off by the beam remnant is (1 — x\ — £2) in DPS and (1 — x) 
in SPS. The results in Fig. [2] show that this remnant momentum is not too different in 
DPS and SPS. Thus, the use of Eq. ([2]) in calculations of event rates at the LHC appears 
adequate as a good first approximation. While available Tevatron data on double parton 



scattering 2J, |25[ are insensitive to possible correlations in x, the greater dynamic range at 
the LHC may make it possible to observe them. 1 

Assuming next that the two subprocesses A(i j —> a b) and B(k I —> c d) are dynamically 
uncorrelated, we express the double parton scattering differential cross section as: 

da DPS = ^- J H ?( x ^ x 2, Va, Mb)#|Vi, Ma, Hb) (4) 

i,j,k,l 

xda i: j(xi, x[, fiA)dd'^ l (x2, x' 2 , iiB)dx\dx2dx\dx ! 2 . 

The symmetry factor m is 1 if the two hard-scattering subprocesses are identical and is 2 
otherwise. In the denominator, there is a factor a c s with the dimensions of a cross section. 
Given that one hard-scatter has taken place, a e s measures the size of the partonic core in 



1 As emphasized in Refs. lg, |22j, even if the approximation in Eq. (2J holds at one hard scale, evolution 



of the parton densities with [i will induce violations at larger scales. 
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FIG. 2: Values of the parton longitudinal momentum fractions x in the DPS and SPS events. Most 
DPS events have low x values. The events used for this plot include the requirements nj et = 4, 
^btag = 2, and the threshold cuts discussed in Sec. II. 

which the flux of accompanying short- distance partons is confined. It should be at most 
proportional to the transverse size of a proton. For the first process of interest in this paper, 
pp — > bbjij 2 , Eq. P| reduces to 

da^ipp -> bb jlh X) = d^PP-bb^da^pp^j^X) ^ (g) 



0"cff 



Tevatron collider data 



24 . |25| yield values in the range a c s ~ 12 mb. We use this value 
for the estimates we make, but we emphasize that the goal should be to make an empirical 
determination of its value at LHC energies. 

In Sec. Ull we present our calculation of the double parton and the single parton contri- 
butions to p p — > b b jx J2 X. We identify variables that discriminate the two contributions 
quite well. In Sec. IHIf we treat the double parton and the single parton contributions to 
4 jet production, again finding that good separation is possible despite the combinatorial 
uncertainty in the pairing of jets. We show in both cases that the double parton contribution 
falls off significantly more rapidly with p£, the transverse momentum of the leading jet. For 
the value of a e g ~ 12 mb and the cuts that we use, we find that, in the region in which it 
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is most identifiable, double parton scattering is dominant for p^! < 30 GeV in b b j\ j 2 at 
LHC energies, and p^ 1 < 50 GeV in 4 jet production. Our conclusions are found in Sec. [TV] 

II. HEAVY QUARK PAIR AND JET PAIR PRODUCTION IN QCD. 

In this section, we describe the calculation of the DPS and SPS event rates for bbjj 
production at the LHC. For our purposes, light jets (denoted by j) are assumed to originate 
only from gluons or one of the four lighter quarks (u, d, s or c) and, as stated above, we 
perform all calculations for the LHC with a center-of-mass energy of = 10 TeV. Event 
rates are quoted for 10 pb _1 of data. 

A. Outline of the method 

The prediction for the DPS event rate is based on the assumption that the two partonic 
interactions which produce the bb and jj systems occur independently (as expressed in 
Eq. (j4])). At leading order, the only contribution is: 

(ij -> bb) <g> (kl -> jj) (6) 

where the symbol <S> denotes the combination of one event each from the bb and the jj 
final states. In an attempt to model some of the effects expected from initial- and final-state 
radiation, we also account for the possibility of an additional jet which is undetected because 
it is either too soft or outside of the accepted rapidity range. Thus, we include several other 
contributions to the DPS event: 

bb{j)®jj , bbj®{j)j , bbj®j{j) (7) 
bb®(j)jj , bb®]{j)j , bb®jj{j), (8) 

where the parentheses surrounding a jet indicate that it is undetected. We compute processes 
such as jj(j) and bb(j) at LO as 3 parton final-state processes. 

The 2 to 3 parton amplitudes for bb(j) [and jj(j)) diverge as the undetected jet (j) 
becomes soft or collinear to one of the other final state partons or to an initial parton. 
The divergences are removed in a full next-to-leading order (NLO) treatment, in which 
real emission and virtual (loop) contributions are incorporated, and the finite bb, bb(j), 
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and bbj contributions are present with proper relative normalization. In the LO parton 
level simulations done in this paper, we employ a cut at the generator level to remove the 
divergences. All the final state objects in the processes listed above are required to have 
transverse momentum > 20 GeV. In this fashion, we model some aspects of the expected 
momentum imbalance between the b and b arising from the 2 to 3 process ij — > bbj, but 
we cannot claim to include the relative normalization between the bb and bbj contributions 
that would result from a full NLO treatment. We leave a complete NLO analysis for future 
work. 

The SPS cross section is computed according to Eq. It receives contributions at 
lowest order from the 2 parton to 4 jet final state process: 

ij -> bbjj , (9) 

and, in the case where a jet is undetected, from the 5-jet final states (computed at LO): 

bb(j)jj , bbj(j)j , bbjj(j). (10) 

We also investigate the possibility of jjjj and jjjj(j) final state contributions to the SPS 
cross section where two of the jets "fake" b jets. We find that the effects from these final 
states are subdominant compared to the processes listed in Eqs. (J9j) and ffTOl) . 

In our numerical analysis, we use the leading-order CTEQ6L1 parton distribution func- 



tions (PDFs) [28| to compute both DPS and SPS cross sections, and we evaluate all cross 
sections using one-loop evolution of a s (/j,). For the renormalization and factorization scales, 
we choose the dynamic scale: 

V 2 = J2 p h + m i> ( n ) 

i 

where px,i is the transverse momentum of the i th jet and rrii = (m-i = 4.7 GeV) for light 
(bottom) jets. In the case of roughly equal values of the transverse momenta pr,i, Eq. ffTTj) 
yields /i ~ 2p T in SPS and /i ~ V2p T in DPS. At LO there is no obviously "right" hard 
scale, and the choice in Eq. ffTTj) seems as good as any other. 

The DPS events are generated as two separate sets of events with Mad- 



graph/Madevent [29] and then combined as described above. For example, at leading order, 
we generate events separately for pp — > bbX and pp — > jjX, and these events are then 
combined as indicated in Eq. (J6j). To increase the speed of the simulations, the SPS events 
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are generated with Alpgen 30] since the SPS processes of interest are hard-coded in Alpgen, 
which contains more compact expressions for the squared-matrix-elements than Madgraph. 

The events accepted after generation are required to have 4 jets rij et = 4 with 2 of these 
tagged as 6's n btag = 2. At the generator level, all the final state objects in the processes listed 
in Eq. ([6]) through Eq. fllOp must have transverse momentum > 20 GeV, as mentioned 
above. Furthermore, at the analysis level, all events (DPS and SPS) are required to pass 
the following acceptance cuts: 

Ptj > 25 GeV, |^| < 2.5 (12) 
p T ,b > 25 GeV, \r] b \ < 2.5 (13) 
ARjj > 0.4, AR bb > 0.4 (14) 

where rji is the jet's pseudorapidity, and AR^ is the separation in the azimuthal angle (0) - 
pseudorapidity plane between jets i and j: 



AR 



»fe) 2 + - &) 2 ■ ( 15 ) 



We model detector resolution effects by smearing the final state energy according to: 

5E a 



(16) 



E ^E/GeV 

where we take a = 50% and b = 3% for jets. To account for b jet tagging efficiencies, we 
assume a 6-tagging rate of 60% for 6-quarks with pt > 20 GeV and \rj b \ < 2.0. We also 
apply a mistagging rate for charm-quarks as: 

e c ^ b = 10% for p T (c) > 50 GeV (17) 

while the mistagging rate for a light quark is: 

e u ,d,s, 9 ^b = 0.67% for p T (j) < 100 GeV (18) 
e uAs>g ^ b = 2% for p T {j) > 250 GeV. (19) 



Over t 
above 



re range 100 GeV < pr{j) < 250 GeV, we linearly interpolate the fake rates given 
31]- 
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B. Properties of SPS and DPS in b b j j 

Having detailed the calculation of the bbjj event rates from DPS and SPS, we now 
discuss some of the distinguishing characteristics of the two contributions. First, however, 
it is important to check that our simulations of DPS events are not introducing an artificial 
correlation between the bb and jj final states. We do this by inspecting the angle $ between 
the plane defined by the bb system and the plane defined by the jj system. If the two 
scattering processes ij — > bb and kl — > jj which produce the DPS final state are truly 
independent, one would expect to see a flat distribution in the angle $. By contrast, many 
diagrams, including some with non-trivial spin correlations, contribute to the 2 parton to 4 
parton final state in SPS, and naively one would expect some correlation between the two 
planes. To avoid possible effects from boosting to the lab frame, we define the two planes 
in the partonic center-of-mass frame. 

We specify the planes by using the three-momenta of the outgoing jets. Then, the angle 
between the two planes defined by the jj and bb systems is: 

cos® = n 3 (j 1 ,j 2 )-h 3 (b 1 ,b 2 ), (20) 

where n 3 (x, y) is the unit three- vector normal to the plane defined by the x — y system. 

The normal is undefined when j\ and ji are back-to-back or b\ and 62 are back-to-back, 
as occurs in a large fraction of the DPS events. Therefore, when cos(p^ x ^ < —0.9, we use a 
different procedure. We use the three-momentum of one of the incoming partons along with 
the three-momentum of one of the outgoing b quarks to define the bb plane. Let qb be the 
three-momentum of an incoming parton, and pb be the three-momentum of the final-state 
b (or b) quark. We then define <p Pb (tb to be the azimuthal angle of the three-vector normal 
to the qb — Pb plane. Note that we use <fi here since the normal to any three- vector and the 
beam-line will be transverse to the beam-line (not the case in the SPS process). In this way, 
the jet which is not used to define the plane is guaranteed to lie in the plane. The plane for 
the jj system is defined in an analogous manner. Finally, the angle between the planes is 
then: 

$ = \<f> Pj , qj ~ Km\ ■ ( 21 ) 

In Fig. El we display the number of events as a function of the angle between the two 
planes. There is an evident correlation between the two planes in SPS, while the distribution 
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FIG. 3: Event rate as a function of the angle between the two planes defined by the bb and jj 
systems. In SPS events, there is a correlation among the planes which is absent for DPS events. 

is flat in DPS, indicative that the two planes are uncorrelated. 

Another interesting difference between DPS and SPS is the behavior of event rates as a 
function of transverse momentum. As an example of this, in Fig. HJ we show the transverse 
momentum distribution for the leading jet (either a b or light j) for both DPS and SPS. 
Several characteristics are evident. First, SPS produces a relatively hard spectrum, and for 
the value of a e g and the cuts that we use, we see that SPS tends to dominate over the full 
range of transverse momentum considered. On the other hand, DPS produces a much softer 
spectrum which (up to issues of normalization in the form of a e s) can dominate at small 
values of transverse momentum. The cross-over between the two contributions to the total 
event rate is ~ 30 GeV for the acceptance cuts considered here. A smaller (larger) value of 
o" c ft would move the cross-over to a larger (smaller) value of the transverse momentum p!^ 
of the leading jet. 

C. Distinguishing variables 

We turn next to the search for variables that may allow for a clear separation of the 
DPS and SPS contributions. Since the topology of the DPS events includes two 2-^2 
hard scattering events, the two pairs of jet objects are roughly back-to-back. We expect 
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FIG. 4: The transverse momentum px distribution of the leading jet in jjbb after minimal cuts. 

the azimuthal angle between the pairs of jets corresponding to each hard scattering event 
to be strongly peaked near A<pjj ~ ~ tt. Real radiation of an additional jet, where the 
extra jet is missed because it fails the threshold or acceptance cuts, allows smaller values of 
The relevant distribution is shown for light jets (non 6-tagged) in Fig. [5^. There is a 
clear peak near A(f)jj = tt for DPS events, while the events are more broadly distributed in 
SPS events. The secondary peak near small Afyj arises from gluon splitting which typically 
produces nearly collinear jets. The suppression at still lower Acpjj comes from the isolation 
cut AR j:j > 0.4. 

The separation of DPS events from SPS events becomes more pronounced if information 
is used from both the bb and jj systems. As an example, we consider the distribution built 
from a combination of the azimuthal angle separations of both jj and bb pairs, using a 
variable adopted from Ref. (25 ]: 

S<p = ^AWubtf + AMjuh)*. (22) 

In Fig. [5]d, we present a distribution in for both DPS and SPS+DPS events. Again, 
as in the case of the A(p distribution, we see that the SPS events are broadly distributed 
across the allowed range of S^. However, the combined information from both the bb and 
jj systems shows that the DPS events produce a sharp and substantial peak near ~ tt 
which is well-separated from the total sample. 

The narrow peaks near A<pjj = tt in Fig. [5^ and near = 1 in Fig. [5Jd will be smeared 
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FIG. 5: (a) The difference A(ft in the azimuthal angles of light jet pairs for DPS and both SPS+DPS 
events. The dijet pairs are back-to-back in DPS events, (b) The variable for DPS and SPS+DPS 
events provides a stronger separation of the underlying DPS events from the total sample when 
compared to Acp for any pair. 



somewhat once soft QCD radiation and other higher-order terms are included in the calcu- 
lation. 

Another possibility for discerning DPS is the use of the total transverse momentum of 
both the bb and jj systems. At lowest order for a 2 — > 2 process, the vector sum of the 
transverse momenta of the final state pair vanishes. In reality, radiation and momentum 
mismeasurement smear the expected peak near zero. Nevertheless, we still expect DPS 
events to show a distribution in the transverse momenta of the jet pairs that is reasonably 
well-balanced. To encapsulate this expectation for both light jet pairs and 6-tagged pairs, 
we use the variable [251 ] : 

q , J_ If \PT(bM\ V ( br(ji,J 2 )| V , . 

PT v^V VIMMI + IMW \\PT(jx)\ + \p T (j2)\) ' { ) 

Here Pr{bi, ^2) is the vector sum of the transverse momenta of the two final state b jets, and 
Pt{jii 32) is the vector sum of the transverse momenta of the two (non b) jets. 

The distribution in S' is shown in Fig. O As expected, we observe that the DPS events 
are peaked near S' ~ and are well-separated from the total sample. The SPS events, on 
the other hand, tend to be far from a back-to-back configuration and, in fact, are peaked 
near S' pT ~ 1. This behavior of the SPS events is presumably related to the fact that a large 
number of the bb or jj pairs arise from gluon splitting which yields a large px imbalance 
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FIG. 6: Distribution of events in SL for the DPS and SPS samples. Due to the back-to-back 
nature of the 2 — * 2 events in DPS scattering, the transverse momenta of the jet pair and of the 
6-tagged jet pair are small, resulting in a small value of S' pT . In (a) we show the S' pT distribution 
for our standard cuts, and in (b) we increase the cut on the transverse momentum of the leading 
jet, p£ > 40 GeV. The fraction of DPS events in the whole sample decreases with increasing pjt. 

In this subsection, we find that extraction of the DPS "signal" for bbjj production from 
the SPS "background" can be enhanced by combining information from both bb and jj 
systems. Our simulations suggest that the variable S' PT may be a more effective discriminator 
than S^. However, given the leading order nature of our calculations and the absence of 
smearing associated with initial state soft radiation, this picture may change and a variable 
such as Sfj, (or some other variable) may become a clearer signal of DPS at the LHC. 
Realistically, it would be valuable to study both distributions once LHC data are available 
in order to determine which is more instructive. In the following, we use the clear separation 
shown in Fig. [6] in our exploration of the distinct properties of DPS and SPS events. 

D. Two-dimensional distributions 

The evidence in Fig. and Fig. [H] for distinct regions of DPS dominance prompts the 
search for greater discrimination in a plane represented by a two dimensional distribution of 
one variable against another. We examined scatter plots involving the inter-plane angle $, 
the jet-jet azimuthal angle difference Acpjj, S</>, and S' . Strong kinematic correlations are 
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evident in the plot of vs. S pT at the level of our leading order calculation, and we observe 
no additional separation of DPS and SPS beyond that evident in Figs. [5] and [6j Likewise, 
there are strong correlations between Acfijj and S^. 

One scatter plot with interesting features is displayed in Fig. [7J The DPS events are 
seen to be clustered near S' pT = and are uniformly distributed in $. The SPS events peak 
toward S' = 1 and show a roughly sin $ character. While already evident in Figs. [3] and [6], 
these two features are more apparent in the scatter plot Fig. [7J Moreover, the scatter plot 
shows a valley of relatively low density between S' ~ 0.1 and ~ 0.4. In an experimental 
one-dimensional $ distribution such as Fig. [3J one would see the sum of the DPS and SPS 
contributions. If structure is seen in data similar to that shown in the scatter plot Fig. [7J, 
one could make a cut at S' pT < 0.1 or 0.2 and verify whether the experimental distribution 
in $ is flat as expected for DPS events. 




FIG. 7: Two-dimensional distribution of events in the variables $ and S' pT for the DPS and SPS 
samples. 



In Fig. HI we show that DPS produces a softer transverse momentum distribution for the 
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leading jet (either a b or light j). In data one would see only the sum of the DPS and SPS 
components in a plot like Fig. LH A scatter plot of S' PT vs. the transverse momentum of the 
leading jet motivates an empirical separation of the two components. In Figs. 0(a) and^b) 
we compare the S' distributions for two different selections on the transverse momentum 

of the leading jet in the bbjj sample. This comparison of the distributions confirms that 
events in the DPS region, defined empirically by the region S' < 0.1 or 0.2, fall off more 
steeply with pl^ than the rest of the sample. It will be important and interesting to see 
whether the selection S pT < 0.1 or 0.2 in LHC data also produces events that show a more 
rapid decrease with pt^ . 

The leading-jet transverse momentum distributions are shown in Figs. [H(a) and EJ^b) 
for two different cuts on S' pT . In both cases, we see that the SPS sample has a broader 
distribution in p^ and that the DPS sample dominates for small enough values of pip . For 
our chosen value of cr e g ~ 12 mb, and for cuts we employ, the crossover points are roughly 
80 GeV for SL < 0.2 and 40 GeV for SL < 0.4. 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 

Leading jet Py (GeV) Leading jet Py (GeV) 

FIG. 8: The distribution in the transverse momentum of the leading jet for (a) S' < 0.2 and 
(b) S' PT < 0.4. As the signal region becomes more dominated by SPS events (i.e. moving from (a) 
to (b)), the resulting distribution becomes harder and shifts the SPS-DPS cross-over from ~ 80 
GeV to ~ 40 GeV. 
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III. FOUR JET PRODUCTION 



In addition to bbjj, we can also ask how important DPS can be for a generic 4j final-state, 
where none of the jets are 6-tagged. In this section, we describe our calculation of the double 
parton scattering and the single parton scattering contributions to the production of a 4j 
final state, for which the cross section is larger. Our exposition can be brief since we repeat 
the procedure described in some detail in Sec. [Til 

A. Outline of the method 

The DPS process for 4j production is topologically equivalent to bbjj. However, in the 
4j system, we lose the 6-tagging ability that reduces the combinatorial background in bbjj, 
and the prospects for isolating and measuring DPS over the SPS background may appear 
less promising. Fortunately, in going from the bb subprocess to the jj subprocess, a much 
larger DPS rate is possible due to the much larger cross section for jj production. As we 
show below, we find that the DPS signature can be extracted in this 4j mode as well. 

The DPS cross section for Aj production receives contributions from the following sub- 
processes at the lowest order: 

j]®]j , bb®jj, (24) 

where both 6-quarks fail the 6-tag. We do not include the bb®bb process due to its relatively 
small rate (~ 0.14 nb). This rate is further reduced by requiring no 6-tags, yielding roughly 
40 events in the 10 pb _1 of luminosity assumed here. 

Following Sec. [TTJ we account for the possibility of an additional jet which is undetected 
because it is too soft or outside of the accepted rapidity range. Thus, we include several 
other contributions to the DPS cross section: 

Jjj ® , (25) 
bbj®j(j) , bb(j)®jj, (26) 
bb®j(j)j , b(b)®jjj , (b)b®jjj. (27) 

where the parentheses surrounding a jet signify that it is not detected. 

The SPS cross section receives contributions at lowest order from the final state: 

jjjj ,bbjj, (28) 
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where both 6-quarks fail the 6-tag, and, in the case where a jet is not detected, from the 
final states: 

bb(j)jj , (b)bjjj , b(b)jjj , (j)jjjj . (29) 

We refer to Sec. [Til for the specification of acceptance cuts and detector resolution, and for 
our treatment of the potential divergences present in the amplitudes for the processes in 
Eqs. (|23)-P3P. 

B. Results 

Similar to the bbjj process, the leading jet in the 4j DPS sample is typically softer than 
in the SPS channels (see Fig. [9]). In this case, again using a e s = 12 mb, we find that the 
cross-over between DPS and SPS dominance occurs near pt — 50 GeV, higher than in the 
bbjj case shown in Fig. HI Improvement in the separation between DPS and SPS in the 4j 

1x1 o 7 
1x1 o 6 

c 

5 
« 

= 1x10 5 

S l 

z 

1x1 4 
1x1 3 

20 40 60 80 100 120 140 

Leading jet P T (GeV) 




FIG. 9: As in Fig. [H but for 4j events. Similar to the bbjj sample, the SPS sample exhibits a 
harder pt spectrum. 

case can be achieved with an analogous version of the S' variable introduced in Eq. (|23l) : 

q > J_ f \Pr(jaJb)\ \ 2 ( \PT(jc,jd)\ \ 2 r „ n x 
PT V2\j \\Pr(j a )\ + \ PT (j b )\J + \\PtUc)\ + \p T {j d )\) ' 1 ' 

Here PT{ja,jb) is the vector sum of the transverse momenta of two final state jets, a and b, 
chosen among the four. The remaining c and d jets are then fixed. This choice is unique if a 
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separation of the two hard interactions is possible. In the bbjj system, the separation into the 
bb and jj subsystems via 6-tagging removed most of the degeneracy (some degeneracy still 
remained via tagging efficiencies or light jet mistagging). In the 4j system, the degeneracy 
can at first glance be problematic as there are 3 possible pairings of the four jets. 




Democratic S D ' 

FIG. 10: The democratic S pT distribution for Aj events shows much more combinatorial background 
than in the bbjj events. Even after accepting two mis-matched jet pairs, we see that the DPS and 
SPS samples can still be separated well. 

One might be tempted to take the pairing of jets which minimizes the value of S' pT . Unfor- 
tunately, this choice places a bias on the distribution that makes it potentially problematic 
to trust the discrimination. Instead, to construct S' we take all three combinations of pair- 
ings, which includes one "correct" pairing and two incorrect pairings in the DPS process. 
This "democratic" S' distribution is shown in Fig. [10] and is re- weighted by 1/3 for proper 
normalization. As in the bbjj case, we see that the DPS distribution peaks near S' PT ~ 0, 
indicative that two back-to-back hard interactions are present. In addition to this expected 
feature, we also see a continuum that extends above S' ~ 0.1, associated with the wrong 
combination taken in the democratic approach. In Fig. [10] we see that DPS produces a sec- 
ondary peak at S' PT ~ 1, not present in the bbjj case in Fig. [6] It appears to arise from the 
wrong pairings of jets associated with the combinatorial background. In these instances, the 
wrong combination of two jets that are close together in AR, meaning that their momenta 
are aligned, can maximize the value of S' . Overall, we see that the DPS peak near S' = 
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provides a good means to separate DPS events from SPS events. 




FIG. 11: As in Fig. El but for 4j events with (a) democratic S' pT < 0.2 and (b) democratic 
S' PT < 0.4. As in bbjj events, as one increases the cut on S' pT , the SPS fraction increases and the 
total distribution is harder. 

As in the bbjj case, we inspect the distribution in the px of the leading jet after cuts 
on the S' Pt variable. Since there are three jet pairings per event, we now require that at 
least one of the three pairings has S' pT in the given window. Due to this softer constraint, 
the hardening of the px spectrum of the leading jet is less dramatic than in the bbjj case 
(e.g. compare Figs. [51 and [TT]) . The crossover of the SPS and DPS contributions occurs near 
80 GeV for S' PT < 0.2 and near 50 GeV for S' PT < 0.4 

IV. DISCUSSION AND CONCLUSIONS 

Our goal is to develop a method to search for a double parton scattering contribution 
in the bbjj and 4 jet final states at LHC energies and to measure the magnitude of its 
contribution relative to the single parton contribution to the same final states. Based on our 
parton level simulations, we find that variables such as S' and that take into account 
information from the entire final state, thereby including both of the hard subprocesses 
in DPS, are more effective at discrimination than variables such as A<f)jj that reflect only 
a subset of the final-state. The enhancement at low values of S' pT shown in Figs. [fU [7] 
and [10] provides a good signature for the presence of double parton scattering. We urge 
experimenters to search for such a concentration of events in data at the LHC. Having 
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found this enhancement, we then suggest that the magnitude of this peak be examined as a 
function of the transverse momentum p^! of the leading jet in the event sample. The double 
parton scattering contribution in the peak region should fall off more rapidly with p^! than 
the rest of the sample. The distribution of events in the region of small values of S' should 
also be examined as a function of the inter-plane angle $ to see whether the flat behavior 
is seen, as expected for two independent production processes. Once these characteristics 
of double parton scattering are established, the data can be used to determine the effective 
normalization a e g, defined and discussed in the Introduction. It will be interesting to see 
whether the values extracted for cr c g are about the same in the b b j j and 4 jet final states 
and how they compare with values measured at the Fermilab Tevatron. 

Once double parton scattering is established in data, and a e g is determined, in a relatively 
clean process such as bbjj, double parton contributions to a wide range of other processes 
can be computed with more certainty about their expected rates at LHC energies. To be 
sure, given the approximations described in the Introduction, some variation in the values 
of (Jeff might be expected and appropriate for different processes and in different kinematic 
regions. The connection of cr e ff with the effective size of the hard-scattering core of the 
proton may mean that a e g will have different values for gg, qq, and qq scattering. 

There are several avenues for future work. Of great importance is the proper inclusion of 



next-to-leading order contributions 



321 ] . They are needed to make more robust predictions 



of the relative normalization of the DPS and SPS contributions, of the shape of the pt 



distribution of the leading jet, and for proper softening of the sharp peaks seen near S' 



PT 



in Figs. [6] and [TUJ and near S<f, = n in Fig. [5b. 

It will also be important to develop joint probabilities H %,k (xi, x 2 , Ha, Hb) that are more 
sophisticated theoretically than the first approximation represented by Eq. (J2J) in which 
parton-parton correlations are absent. A valuable development in this direction are the 
studies presented in Refs. |16l. l22j]. 

Double parton contributions are potentially relevant for a wide range of standard model 
many alread y co nsidered in the literature 



processes 
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131, 



251 ] , a nd they may also feed pertinent standard model 



backgrounds to new physics processes [33j. They could be an issue in studies of Higgs boson 
production in weak-boson-fusion since the "forward" jets could come from a second hard 
subprocess. 
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